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The hadronic decays of Xci — >■ VV and Xc2 — > VP are supposed to be suppressed by the helicity 
selection rule in the pQCD framework. With an effective Lagrangian method, we show that the 
intermediate charmed meson loops can provide a mechanism for the evasion of the helicity selection 
rule, and result in sizeable decay branching ratios in some of those channels. The theoretical 
predictions can be examined by the forthcoming BES-III data in the near future. 
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I. INTRODUCTION 

_ ^ The exclusive decays of heavy quarkonium have been an important platform for studying the nature of strong 
. interactions in the literature [JQ since the discovery of quantum chromodynamics (QCD). In this energy region 
relatively large energy scale (~ TOc,Wb) is involved such that some perturbative QCD(pQCD) asymptotic behaviors 
^ , can be expected, for instance, the so-called helicity selection rule By studying the manifestation of the helicity 
selection rule in heavy quarkonium decays, we expect to gain insights into the property of QCD in the interplay of 
\ the perturbative and non-perturbative regime. 

■ We briefly review this powerful method that is elaborated on in Ref. [2]. For a charmonium meson Jcc decaying 
into two light mesons hi and /i2, the perturbative method gives the asymptotic behavior of the branching ratio as 
I— I . follows 

^ '. 

where A, Ai, and A2 are the helicities of the corresponding mesons. This is the result of the pQCD method to leading- 
twist accuracy; i.e. only the valence Fock state (here it is cc) is considered. It is obvious that the leading contribution 
I comes from when Ai + A2 = 0, (while the higher twist will be suppressed by at least a factor of Aqp^/ttt,^.) while the 
helicity configurations which do not satisfy this relation will be suppressed. 

An alternative description of this selection rule is to depict it with the quantum number "naturalness" , which is 
ly-j , defined as cr = P(— 1)'^, where P and J are the parity and spin of the particle, respectively. The selection rule then 



(N 



requires that 



^initial = criCTz • (2) 

That is, the initial-state naturalness should be equal to the product of those of the final states. We can comprehend 
• • this in such a way that, if cr™'*'"' o'i<J2, to keep the parity conservation and Lorentz invariance, the amplitude will 
. 5^ contain a Levi-Civita tensor e^,yap , which will be contracted with the polarization vectors and momenta of the involved 
. mesons. For instance, for the process Xci W where VV represent a pair of vector mesons, in the rest frame of 
' Xci J the non- vanishing covariant amplitude will require one vector meson be transversely polarized and the other one 
5^ I is longitudinally polarized. This leads to Ai -f A2 7^ 0, which will violate the helicity selection rule and is supposed to 
be suppressed. Another well-known example is the process J/ijj — > VP, where the non-vanishing covariant amplitude 
also violates the selection rule. 

The helicity selection rule has been used for studying some exclusive decay processes of heavy quarkonia. In 
Refs. [3, [1, 0|i some allowed or suppressed processes are explicitly listed. Interestingly, with the accumulation of the 
experimental data, more and more observations suggest significant discrepancies between the data and the selection- 
rule expectations. For instance, the decays of J/ip — > VP and rjc — VV would be suppressed by this rule. In reality. 
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they are rather important decay channels for J/ip and ?7c, respectively ]9^. One possible reason why the perturbative 
method fails here could be that although the mass of the charm quark is heavy, it is, however, not as heavy as pQCD 
demands. Therefore, it is not safe to apply the helicity selection rule to charmonium decays, while the situation may 
be improved in bottomonium decays. Taking into account these issues, it implies that there might be some other 
mechanisms that can contribute to those helicity-selection-rule-forbidden processes, such as higher order contributions, 
final state interactions, or some other long-distance effects. Theoretical discussions on these mechanisms in B-meson 
or heavy quar konium decays can be found in the literature 

In Ref . [iSl il6j , the role played by the intermediate charmed meson loops in the understanding of the so-called 
"p-TT puzzle" and ■0(3770) non-DD decays into VP was studied. It shows that the long-distance interactions via the 
intermediate meson loops provide a possible mechanism for understanding the deviations from pQCD expectations. 
In the present work, we will study the P-wave charmonium decays, i.e. Xci ~^ and Xc2 VP, which are 
suppressed by the helicity selection rule but may be enhanced by the intermediate charmed meson loop transitions. 
The branching ratio of Xci ^*o^*o order of 10~^ has been measured by experiment ^8|]. According to the SU(3) 
flavor symmetry, it is justified to predict some other channels as pp, ljlo and 00, of which the branching ratios may 
also be sizeable. Therefore, it is interesting to investigate the intermediate meson loop transitions as a mechanism 
for evading the helicity selection rule. We also note that the process Xc2 — > VP will be further suppressed by the 
approximate G-parity or isospin ([/-spin for strange mesons) conservation, and the Xc2 decays into neutral VP with 
fixed C parities are forbidden by C-parity conservation. By comparing the results with the SU(3) flavor symmetry 
expectation, we anticipate that the non-perturbative mechanism can be highlighted. 

It is worth noting that for the P-wave charmonium exclusive decays, the next higher Fock state ccg^ i.e. the so- 
called color-octet state, will contribute at the same order as the color-singlet state cc in the framework of perturbative 
method 0, [13423 • This scenario may share the same intrinsic physics with the intermediate meson loop transitions, 
but based on different view angles. Namely, one description is at the quark-gluon level, while the other is at the 
hadron level. Taking into account the relatively low energy scale, it is not easy to handle these processes with the 
factorization method for the quark-gluon interactions. Some qualitative discussions can be found in [ll|, [I^, [3| and 
references therein. 

This paper is arranged as follows; In Sec. II, the formulae for the intermediate meson loops as a long-distance effect 
in the charmonium decays will be given by an effective Lagrangian method. The numerical results and discussions 
will be given in Sec. III. The conclusion will be drawn in Sec. IV. 



II. LONG-DISTANCE CONTRIBUTION VIA INTERMEDIATE MESON LOOPS 

In Figs. [U [5] and 131 the intermediate meson loops are plotted for xa VV, Xc2 PT^ and Xc2 — >■ K*K, respectively. 
The relevant effective Lagrangians that will be used are based on the heavy quark symmetry and SU (3) symmetry 
[isl [l^ [2^ . The spin multiplet for these four P-wave charmonium states are expressed as 

P& = (^) {x'i2l^ + ^e.^^.pv'^^^x:, + ^(7" - «'')XcO + Kl,^ (^i^) . (3) 
The charmed and anti-charmed meson triplet read 

Hu = [2?:% -25.75], (4) 

H2^ = [A*^M-^.75](^), (5) 

where T) and V* denote the pseudoscalar and vector charmed meson fields respectively, i.e. P*-*-* — 
(D°(*), £>+(*), D^^*^). Then the effective La grangian that describes the interaction between the P-wave charmonium 
and the charmed mesons reads 

£1 = igiTr[P,^sH2^j^Hu] + H.c. (6) 
The effective Lagrangians describe the couplings of charmed mesons to light hadrons read 

Cvvv = -igDDv'D,'^^,VJ{V^'),J, (7) 
Cvvv = -2fD^DVt^uc.p{d''V''),j{V,d Vf-Vfd V,), (8) 
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FIG. 1: Triangle loop diagrams that describe the long-distance contributions in Xci ^ V1V2, where V1V2 represent K* K* , pp, 
ujuj, and 00, respectively. 



1 ~* 1/ *f3 

CvvV ~ ■^gD*D'-p£p,ual3'Pi'^d''Vij d T^* , 

where V and V are the pseudoscalar octet and vector nonet, and we take the following conventions 



p 

K* 



(9) 
(10) 

(11) 



(12) 



K* 



V2 



K- 



h 



(13) 



Xci ^ vv 



Figure [T] illustrates the long-distance contributions for Xci — >■ VV . One notices that there is no such a vertex of 
XciD* D* in Fig. [TJ It is because to conserve parity and keep Lorentz invariance, the relative angular momentum 
between D*D* should be L = 2, and the total spin 5 = 2. But such kind of coupling is not present in the 
expansion of the effective Lagrangian Ci. For Fig. [Tfa) and (b), we take the convention of the momenta as Xciip) ~^ 
D{qi)D* {q2)[D'^*\q)] — Vi(pi)V2(p2), where the exchanged particle between D and D* is indicated in the square 
bracket. Figure [TJc) and (d) will give the same contribution as Fig. [IJa) and (b). For simplicity, we just write down 
the amplitude Aiia and AAib and take the final states p~ as an example: 



Mia = '2igDD'x.i9DDvfD'DV<^l''^<^r<^*'' 



92, 



A /3 / 
92 92/" 



M-ib — '2ig DD'' Xci 9 D DV fo'DV^^"^ £i'^£2^ 



DaDiD2 
d^q 



(14) 



(2^) 



X e^<Ta/3Pi (9? + 9") [gD*D'v{q2r ~ 9T)g75 + D* D'V {p2SgT-i ^ P2',gSr)] 

1 2^ 



X (/^ - qPqyml,){g^' - qUH^l^) x 



(15) 



where Do = q 
obtained similarly. 



m|,, Db = q^ — fn'o*, Di = q\ — rnj-,, and D2 = (/I ~ rnj-,,. The results of other channels can be 
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FIG. 2: Triangle loop diagrams that describe the long-distance contributions in Xc2 — > P^t^ ■ The diagrams for Xc2 — >■ p vr^ 
are implicated. 




FIG. 3: Triangle loop diagrams that describe the long-distance contributions in Xc2 — > K K . The diagrams for Xc2 K K 
are implicated. 



Since the mass of Xci is under the threshold of DD* , intermediate mesons D and D* can not be on-shell simulta- 
neously. We phenomenally introduce a form factor J-(q^) as has been done in Refs. [islll^. 



Hi') = n 



A 



9i 



A? 



(16) 



where qi = q, qi, q^- The cut-off energy is chosen as hi = rrii + aAgcu, and to; is the mass of the corresponding 
exchanged particle. This is somehow different from the form factor adopted in Refs. [H, We should mention 
that the form factor is necessary for killing the divergence of the loop integrals, although it will also give rise to the 
mo del- dependent aspects of the calculations. We will discuss this in detail later. 



B. 



Xc2 



VP 



As mentioned previously, the decay Xc2 — >■ VP suffers not only from the suppression of the helicity selection rule, 
but also from the approximate G-parity or isospin/C/-spin conservation. However, because of the relatively large mass 
difference between the u/d quark and s quark, the intermediate meson loops may still bring in sizeable branching 
ratios for Xc2 — ^ KK* + c.c Next, we will consider, Xc2 p~^iT~ + c.c. and Xc2 — >■ K*K + c.c. The Xc2 decays into 
neutral VP with fixed C-parity, such as p'^n'^ and lot] etc., forbidden by C-parity conservation. The loop diagrams 
for these two processes are presented in Fig. [2] and Fig. [3l respectively. The convention of the momenta follows that 

of Fig. m 

The relative signs between the following amplitudes are opposite, i.e. A^2a and A42c, M.2b and M.2d: A^3a and 
A^3c, and Ai^b and AAm- This leads to destructive interferences between these amplitudes, and is a reflection of the 
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approximate isospin or [/-spin invariance. We explicitly list the amplitudes, A^2a and A^26, 

X {g^^ - qU/ml,){g^^ - q^q^/ml,) ^ ^(g^), (17) 
1 /" 

A^26 = -^igD'D'Xc2gD'D'P€}^{'e'p+ J J^<^paapP2{q" - q2) 

X (/^ - 9?g7M!..)(.9'"' - 9;'92'M!.0(/'' - q'q'/^^l')jr^Hq') , (18) 

while the others can be obtained similarly. As for the polarization sums of e|^^ , we follow the expressions in Ref. [2^ . 

III. NUMERICAL RESULTS 

Before proceeding to the numerical results, we first discuss the parameters, such as the coupling constants, in the 
formulation. In the chiral and heavy quark limit, the following relations can be obtained [isl, [2^: 

(3gv J, fD^D^'V ^gv 

gODV = gO'D'V = —f=-i JD'DV — — —7=-, 

V2 ruD' V2 

rup go'DTT 2 

gv - go'D'TT — . ■ - -^g, 



gO'-D^K — \ gO'DiT, gO'DK = \ , [i-V) 

V mo V "^i?* 

where (3 and A are commonly taken as (3 ~ 0.9, A ~ 0.56 GeV^^ 0, [l3j HBl j while /tt is the pion decay constant. With 
the measured branching ratio of D* — > Dtt by CLEO-c, the coupling g is determined as 5 = 0.59 27]. In the heavy 
quark limit, the expansion of the effective Lagrangian Ci leads to relations for the P-wavc charmonium couplings to 
the charmed mesons as follows: 



goD'x.i 2\^giy/m]ym]pifr^ 



where gi has been related to the XcO decay constant fxao- '^^^ be approximately determined by the QCD sum rule 
method, i.e. /^.o - 0.51 GeV [H]. 

The form factor parameter a generally cannot be determined from the first principle. It is usually taken to be order 
of unity, and depends on the particular process. In this work, since the branching ratio of Xci ^*o^*o j^g^g been 
measured in experiment Q, we will adopt the data to constrain the form factor parameter. 

The loop integrals are calculated with the software package LoopTools [1^ . In Tables U and |TT1 we display the 
numerical results for the branching ratios of Xci — >■ VV and Xc2 VP at a = 0.3 ^ 0.33, which correspond to the 
lower and upper bounds of BR{xci ^ K*°K*°) = (1.6 ± 0.4) x lO"^ ^8j_ 

For Xcl — > W, the branching ratio of the pp channels is significant, while the ujlu and (pcj) channel are relatively 
small. We also include the predictions from the SU(3) flavor symmetry as a comparison. The parametrization is given 
in the Appendix. It is interesting to see that the results from the SU(3) flavor symmetry are basically compatible with 
those given by the intermediate meson loops except that the branching ratio for (fx/) is larger. It should be an indication 
for the SU(3) flavor symmetry breaking. With the SU(3) symmetry breaking parameter R ~ fyr/ /k = 0-838 [s*], the 
branching ratios agree with the intermediate meson loop results pretty well. 

Some insights into the transition mechanisms can be gained here: 

i) The SU(3) flavor symmetry breaking R — 0.838 is consistent with the loop transition behavior. Note that the 
intermediate Ds and D* pair has higher mass threshold, and the production of the 4)(j) will be relatively suppressed 
in comparison with the non-strange pir and lou apart from the final-state phase space differences. This corresponds 
to the flavor symmetry breaking at leading order. 
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BR (xlO~*) 




PP 




H 


Exp. data 


16 ±4 








Meson loop 


12 ~ 20 


26 ~ 54 


8.7 ~ 18 


2.7 ~ 4.6 


SU(3)(i?= 1) 


16.0 


26.8 


8.8 


6.8 


SU(3)(i? = 0.838) 


16.0 


32.0 


10.6 


4.0 



TABLE I: Branching ratios for Xci — ^ W predicted by the intermediate meson loop transitions in the range a = 0.3 ~ 0.33 
corresponding to the measured lower and upper bound of BR{xc2 — > K*''K*°) Two results from the SU(3) flavor symmetry 
relation are presented with R — 1 and R — 0.838. The long-dashed line means that experimental data are unavailable. 



BR(xlO"^) 


K*°K° + c.c. 


A'*+A' + c.c. 


P^TT + C.C. 


Meson loop 


4.0 ~ 6.7 


4.0 ~ 6.7 


(1.2 ~ 2.0) X 10"^ 


Exp. data 









TABLE IL Branching ratios for Xc2 — VP predicted by the intermediate meson loop transitions in the same range of a = 
0.3 ~ 0.33. The long-dashed line means that experimental data are unavailable. 

ii) The consistencies support the idea that the intermediate meson loops provide a nature mechanism for the evasion 
of the helicity selection rule as a long-distance transition. Since the Xc0,i.2 are P-wave states, the short-distance 
transition probes the first derivative of the cc wavefunction at the origin, which however, would be suppressed in 
Xci VV and Xc2 — > VP due to the helicity selection rule. By annihilating the cc at long distance via the 
intermediate meson loops, the helicity selection rule is then evaded. 

iii) This phenomenon is slightly different from the S'-wave charmonium decays, such as J/4' Ptt etc, where it is 
not easy to separate out the long-distance transitions from the short-distance ones arising from a possible anomalous 
component of wavefunction [15|, i30, i31|| • 



For Xc2 — >■ VP, since there are no data available at this moment, our predictions are based on the same form 
factor parameter for Xci VV. As listed in Tab. |lTl the accessible channels are only K*K + c.c. and p'^Tr~ + c.c. 
Interestingly, the branching ratio for K*K + c.c. turns out to be sizeable, and the charged pn is found much smaller 
than the K*K + c.c. channel. Qualitatively, this is because that the cancellations between (a) and (c) (and similarly 
between (b) and (d)) in Figs. [2] and [3] are rather different. Namely, the K*K + c.c. channel experiences large [/-spin 
symmetry breakings due to the significant mass difference between d and s quark. In contrast, the pi: channel is 
originated from the u-d mass difference. This unique result can be examined by BES-III as a test of our model. 



The sensitivities of the calculation results to the form factor parameter a are presented in Figs. Blan d [5 ] One notices 
that the values for a in this work are relatively smaller than those adopted in some other worksfisl [itI [Tsj . This is 
acceptable since in this work the form factor takes care of off-shell effects arising from those three intermediate mesons, 
instead of only the exchanged one in the rescattering [l^ [13, [3 • Given the inevitable mo del- dependence introduced 
by the form factor, what turns to be relatively stable and less model-dependent is the relative branching ratio fractions 
among those decay channels within the adopted range of a. As a consequence, the theoretical predictions for other 
decay channels can be better controlled by the data for Xci K*^K*^ . 

IV. CONCLUSION 

In this work, we have discussed how the long-distance transitions via the intermediate meson loops would contribute 
to the processes xa — ^ VV and Xc2 VP, which are supposed to be suppressed according to the helicity selection 
rule in QCD. With an effective Lagrangian method with heavy quark and chiral symmetry, this helicity-selection-rule 
evading mechanism is quantified. Although there are still relatively large uncertainties arising from the form factor 
parameter, we argue that the fewer sensitivities of the branching ratio fractions among the accessible channels would 
provide a better control of the theory predictions. Reasonable ranges of the predictions are obtained. 

We also compare the results from the intermediate meson loops with the expectations of the SU(3) fiavor symmetry, 
and find that they are in a good agreement with each other. In particular, the SU(3) flavor symmetry breaking would 
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FIG. 4: a-dependence of the calculated branching ratios for Xci — >■ VV. 
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FIG. 5: a-dependence of the calculated branching ratios for Xc2 — > VP. 

lead to a suppression of Xci ^ in comparison with the lolu. It can be well-understood by the heavier mass of the 
intermediate Dg (Dl) than the D {D*) state. Namely, the mass threshold of the DgD* +c.c. is higher than DD* +c.c. 
In Xc2 — > VP, the predicted branching ratio for the K*K + c.c. channel is at order of 10~^, while the charged pir 
channel is small. The suppression on the charged pn (in comparison with the K*K + c.c.) can be comprehended as 
a consequence of the larger effects due to the [/-spin symmetry breaking rather than the isospin symmetry breaking, 
i.e. TTis — rrid >> rud — rriu. 

In brief, we emphasize that the P-wave charmonium decay should be ideal for examining the evading mechanisms 
of the helicity selection rule. Our predictions can be examined by the high-statistics Xcj production in the BES-III 
experiment [y, [s^l ■ 
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Appendix 



We adopt a simple parametrization [33l - l35j to study Xci — >■ VV based on the SU(3) flavor symmetry. By assuming 
that H represents the potential for the cc annihilating into gluons and then hadronizing into {qq)vii<l<l)v2i we define 
the transition amplitude strength as 



90 = {{qq)vi{qq)v^H\x: 



cl) 



(21) 



where q (q) is a non-strange quark (anti-quark). Considering the SU(3) flavor symmetry breaking, we introduce 
parameter R which describes 



R = {iqs)vAsq)v2\H\xci) / {iqq)vi{qq)v2\H\xc 



and 



R 



2 _ 



((ss)vi {ss)v2 \H\xci) / {{qq)vi {qq)v2\H\x^ 



cl) 



(22) 



(23) 



where the exchange of Vi and V2 is implicated, and the SU(3) flavor symmetry is recognized by i? = 1. This parameter 
can be related to the ratio of the tt and K meson decay constant, i.e. 



R — Itt/ Ik , 



which gives i? ~ 0.838 |8|. 

A commonly adopted vertex form factor is also applied. 



J-^(Pi)^|pi|^^exp(-pf/8/3^), 



(24) 



(25) 



where pi is the three- vector momentum of the final-state meson Vi in the rest frame of Xci, and L is the relative 
orbital angular momentum between Vi and V2. As discussed earlier, L = 2 is required by parity conservation and 
Lorentz invariance. We adopt /? = 0.5 GeV, which is the same as in Refs. [1^ I33.-.36.] . The partial decay widths are 
given as follows: 



r 



= r 



Ipil 



p 7r"f 



247rM2 

Acl 



50-^' (Pi) , 



r^.o/f.o — Tx' + K' ^K--K- 

IPlI 

247rM2^^ 



iPil 



247rM2 , 

Acl 



ffo-^'(Pi) 



iPil 



247rM2^^ 



(26) 
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